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Abstract 

Traditional Chinese medicine (TCM) is one of the oldest East Asian medical systems. The present study adopted a systems 
biology-based approach to provide new insights relating to the active constituents and molecular mechanisms underlying 
the effects of dragon's blood (DB) tablets for the treatment of colitis. This study integrated chemical analysis, prediction of 
absorption, distribution, metabolism, and excretion (ADME), and network pharmacology. Firstly, a rapid, reliable, and 
accurate ultra-performance liquid chromatography-electrospray ionization-tandem mass spectrometry method was 
employed to identify 48 components of DB tablets. In silico prediction of the passive absorption of these compounds, 
based on Caco-2 cell permeability, and their P450 metabolism enabled the identification of 22 potentially absorbed 
components and 8 metabolites. Finally, networks were constructed to analyze interactions between these DB components/ 
metabolites absorbed and their putative targets, and between the putative DB targets and known therapeutic targets for 
colitis. This study provided a great opportunity to deepen the understanding of the complex pharmacological mechanisms 
underlying the effects of DB in colitis treatment. 
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Introduction 

Traditional Chinese medicine (TCM) is one of the oldest 
medical systems of health care, and it has been used in East Asian 
countries such as China, Japan, and Korea for thousands of years 
[1,2]. However, acceptance of traditional Chinese medicines 
within W estern biomedical practice has been restricted by a lack of 
knowledge of the active compounds involved, and their therapeu- 
tic mechanisms of action [3] . TCM is characterized by the usage 
of multi-component, multi-target agents that collectively modulate 
molecular networks. Recently, network pharmacology has pro- 
vided a means to improve understanding of the molecular 



mechanisms underlying the therapeutic effects of traditional 
Chinese medicines [3—5]. TCM chemical databases such as 
TCM Database@Taiwan [6] and HerbBioMap database [7] and 
some literature generally provide the main sources of information 
on the chemical profiles of traditional Chinese medicines. 
However, these are not always accurate because the chemical 
profiles of traditional Chinese medicines can vary significantly, 
depending on the geographical origin of the materials used, the 
harvest time, pretreatments, and manufacturing processes em- 
ployed. In addition, most TCM formulations are taken orally, 
estimation of intestinal absorption and cytochrome P450 metab- 
olism provide more in-depth insights into their therapeutic 
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Figure 1. A schematic diagram of this systems biology-based investigation into the molecular mechanisms of DB tablet for colitis 
by integrating chemical analysis, ADME prediction and network pharmacology. 

doi:10.1371/journal.pone.0101432.g001 



mechanisms [8,9]. However, experimental determination using 
these systems can be costly and time-consuming. To obtain a rapid 
estimation of human absorption and metabolism, high throughput 
screening methods using Caco-2 cell monolayers and P450 
enzymes provide the most advanced in vitro approaches to 
assessing new chemical entities [10-12]. Thus, the development of 
rapid high-throughput approaches to studying TCM compositions 
and mechanisms, combining chemical analysis, prediction of 
absorption, distribution, metabolism, and excretion (ADME), and 
network pharmacology is required. 

According to the Pharmacopoeia of the People's Republic of 
China, commercially available longxuejie or dragon's blood (DB) 
includes plant resins from four species; Dracaena spp., Daemonor- 
ops spp., Croton spp., and Plerocarpus spp., and has long been 
used as an ethnomedicine in China to invigorate blood circulation 
in the treatment of traumatic injuries, blood stasis, and pain 
[13,14]. DB has also been used to treat chronic colitis in recent 
decades and many clinical studies have reported good therapeutic 
effects [15-17]. Chemically, flavonoids, phenols, steroides, and 
terpenoids have been identified as the main constituents of DB 
[18-22]. Some analytical methods have been developed to identify 
and characterize DB composition and assess its quality [23-25]. 
Pharmacologic studies of these components have identified a wide 
variety of actions, including anti-inflammatory, antiulcer, antimi- 
crobial, hemostatic, and analgesic activities, which could contrib- 
ute to their efficacy in the treatment of chronic colitis [26-28]. 
However, the precise active constituents of DB that have benefits 
in chronic colitis, and their molecular mechanisms of action, are 
still unclear. In particular, the lack of elucidation of the ingredient- 
target network has hindered understanding of the molecular 
mechanisms of DB in chronic colitis treatment. 

The present study employed high-throughput analysis, in silico 
ADME models, and a network pharmacology technique to 
investigate the active constituents of DB and their actions on 
molecular networks in chronic colitis, as shown in Figure 1. 

Materials and Methods 

1. Reagents and chemicals 

High performance liquid chromatography (HPLC)-grade ace- 
tonitrile, formic acid, and methanol were obtained from Merck 
(Darmstadt, Germany). Water was purified using a Milli-Q_ system 
(Millipore, Billerica, MA, USA). 7,4'-Dihydroxyflavan, loureirin 
A, loureirin B, dracaenin A and pterostilbene standards were 
purchased from the National Institute for the Control of 
Pharmaceutical and Biological Products (Beijing, China). The 
purities of all standards were no less than 98% and suitable for 
liquid chromatography-tandem mass spectrometry (LC-MS/MS) 
analysis. DB enteric-coated tablets were supplied from Yunnan 
Datang Hanfang Pharmaceutical Co. Ltd. (Yunnan, China). 

2. Preparation of samples and standard solution 

Six different batches of DB enteric-coated tablets were 
pulverized with a 60 mesh, respectively. Each sample of this 
powder (0. 1 g) was weighed precisely and ultrasonically extracted 
in 10 ml methanol for 20 min at room temperature. The solution 
was centrifuged at 12000 rpm for 10 min, and then filtered 
through 0.22 urn nylon membrane filters. The filtrate was 
analyzed directly by UPLC-ESI-MS/MS, as described in section 
2.3. At the same time, a stock solution containing five standards 



(7,4'-dihydroxyflavan, loureirin A, loureirin B, pterostilbene, and 
dracaenin A) was prepared in methanol. All solutions were stored 
at 4°C prior to analysis. 

3. Instrument and UPLC-ESI-MS/MS conditions 

UPLC was performed on a Dionex UltiMate 3000 system 
(Dionex Corporation, Sunnyvale, USA) equipped with a quater- 
nary pump, an online vacuum degasser, an autosampler, and an 
automatic thermostatic column oven. A Thermo Hyper Gold C18 
column was used (100x2.1 mm, 1.9 u.m) at 30°C with a flow rate 
of 0.4 ml/ min and an injection volume of 5 u.1. The mobile phase 
was a mixture of 0. 1 % formic acid in water (A) and acetonitrile (B). 
The mobile phase gradient program was 5% B, 0-3 min; 5-95% 
B, 3-20 min; 95% B, 20-25 min. High-resolution accurate-mass 
full scan LC-MS and LC-MS/MS analyses were performed using 
Thermo Q;Exactive (Thermo Fisher Scientific, Bremen, Ger- 
many). Full scans were acquired in the mass analyzer at 100- 
1500 m/ z with a resolution of 70000 in both positive and negative 
ion modes, and MS /MS scans were obtained with a resolution of 
1 7500 using a normalized collision energy of 30% for high-energy 
collisional dissociation fragmentation. Based on the best response 
for most compounds, the final parameters were set as follows: 
spray voltage, 3.2 kV in the positive mode and 3.0 kV in the 
negative mode; capillary temperature, 350°C; sheath gas pressure, 
35 arbitrary units; auxiliary gas pressure, 10 arbitrary units; and 
heater temperature, 300°C. Xcalibur 2.2 software (Thermo Fisher 
Scientific) was used for data acquisition. 

4. Prediction of absorbed constituents and their 
metabolites using in silico ADME models 

Structural information (*.mol or *.sdf files) relating to the DB 
components(identified by UPLC-ESI-MS/MS) were downloaded 
from ChemSpider (http://www.chemspider.com/). ADME eval- 
uation of these constituents was carried out using ACD/Percepta 
software 5.07 (ACD/Labs, Toronto, Canada), including the 
passive intestinal permeability of Caco-2 module and the P450 
regioselectivity module to predict their oral bioavailability. The 
apparent permeability coefficient (Papp) of one constituent was 
greater than 9.0x 10 6 cm/s, which indicated good absorption 
characteristics. The score and the reliability of metabolic reaction 
were set at greater than 0.6 and 0.5, respectively, in order to 
increase prediction credibility. 

5. Known therapeutic targets for colitis treatment 

Known therapeutic targets were obtained from the DrugBank 
database [29] (http://www.drugbank.ca/, version: 3.0). We only 
included drug-target interactions where the drugs were approved 
by the Food and Drug Administration, USA (FDA) for the 
treatment of colitis, with human gene/protein targets. To facilitate 
data analysis, all protein identification codes were converted to a 
common UniProtKB-Swiss-Prot code, and detailed target infor- 
mation is provided in the supplementary Table SI. 

6. PPI data collection 

PPI data were imported from eight existing PPI databases, 
including the Human Annotated and Predicted Protein Interac- 
tion Database (HAPPI) [30], Reactome [31], Online Predicted 
Human Interaction Database (OPHID) [32], InAct [33], Human 
Protein Reference Database (HPRD) [34], Molecular Interaction 
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Figure 2. UPLC-ESI-MS/MS of DB in positive ion mode (A) and negative ion mode (B). 

doi:10.1371/journal.pone.0101432.g002 



Database (MINT) [35], Database of Interacting Proteins (DIP) 
[36], and PDZBase [37]. Detailed information on these PPI 
databases is provided in the supplementary Table S2. 

7. Pharmacological mechanism analysis 

7.1. Prediction of putative DB targets. As described in our 
previous study [38], we used the Drug Similarity Search tool in the 
Therapeutic Targets Database [39] (TTD, http://xin.cz3.nus.edu. 
sg/group/cjttd/ttd.asp, Version 4.3.02, released on Aug 25th 
2011) to screen drugs similar to DB via structural similarity 
comparisons. We only selected drugs with a high structural 
similarity score (>0.85, similar to very similar) with the potentially 
absorbed DB constituents and their metabolites. The therapeutic 
targets of these similar drugs were also included as putative DB 
targets. 

7.2. Network construction and analysis. The absorbed 
DB components and their metabolites, their putative targets, and 
known therapeutic targets for colitis treatment were used to 
construct a chemical component-putative target network, and a 
putative DB targets-known colitis therapeutic targets PPI network, 
respectively. PPI data were obtained from eight existing PPI 
databases, as described in section 2.6. Navigator software (Version 
2.2.1) and Cytoscape (Version 2.8.1) were utilized to visualize the 
networks. 

For each node "i" in these two networks, we defined four 
measures of topology: (1) "degree," defined as the number of links 
to node i; (2) "betweenness," defined as the capacity of node i to be 
located in the shortest communication paths between different 



pairs of nodes in the network [40] . The betweenness centrality of 
node i is computed as following formula: 

U ((N-\)(N -2)/2) 

where a st denotes the number of the shortest paths between node s 
and node t in the PPI network, a st (i) denotes the number of the 
shortest paths across node i between node s and node t, and N is 
the total number of nodes in the PPI network. This property 
correlates more closely with essentiality than connectivity, 
exposing critical nodes that usually belong to the group of scaffold 
proteins or proteins involved in crosstalk between signal pathways 
(called bottlenecks) [41]. (3) "closeness," defined as the inverse of 
the sum of node i distances to all other nodes. Closeness can be 
regarded as a measure of how long it would take to spread 
information from node i to all other nodes sequentially. Degree, 
betweenness and closeness centralities correlate with a protein's 
topological importance in the PPI network [40]. (4) K-core 
analysis is an iterative process in which nodes were removed from 
the network in order of the least-connected [42]. The core of 
maximum order is defined as the main (highest) k-core in the 
network. A k-core sub-network of the original network can be 
generated by recursively deleting vertices from the network whose 
degree is less than k. This results in a series of sub-networks that 
gradually reveal the globally central region of the original network. 
On this basis, "k value" is used to measure the centrality of node i. 



PLOS ONE | www.plosone.org 



4 



July 2014 | Volume 9 | Issue 7 | e101432 



Uncovering the Molecular Mechanisms of DB in Colitis Systemically 



z 

■D 

c 

3 
O 

a 
£ 
o 
u 



-Q 

ru 

T3 
QJ 

fD 
O 



QJ 
T3 



QJ 
Q. 



T3 



re 



I I 2 



£ I 



_c o 

+7^ -S 



O ^ 



q 
x 



° 1 



o o o o 



0 0)01 



Ol Ol Dl 



CO d ^ 



ON <— vO CO 



u-i co m m 



m u-i <— lo i— 



«— co co in 
m rsi (N cm 



cr> ro co t— 



vo u-i on co 



m u-i m lti 



co on rs 
o >- o 
O CO on 



co on co 



O «- t- •- 



r- ON fN 



O >— >— 



on on m m 



m i— m i— 



(N m in rM 



0> vO CO CO 



on d o> 



m m -— 



t— co vo on 



o on m 

(N in q IN 

IS r ffi ui 

<N co i— 



in ir> m in 



<— \D «— 



*D IN vO CO 



m m i— 



VO CO ON 



i— vO rsi 



i— co co in 
rn N (N IN 



m u-i u-i in 

(N IN IN IN 



ON ON in ro 



i — ON IN 



O i— _ ■-; 

m" ro 

is is o o 

m m in in m m 



ON VO co co 



on on on 



O i— ON ro 

<N (N p (N 

is" ro vd 

(N i— CO i— 









■D 






01 






3 












<U 

at 






S 










CO 


i 




ro 




CO 


CO 



rs 


yD 


CO 


fN 




is 


ON 


oS 


oS 



«- i— r* 



r\l r\i r\i rN 



PLOS ONE | www.plosone.org 



5 



July 2014 | Volume 9 | Issue 7 | e101432 



Uncovering the Molecular Mechanisms of DB in Colitis Systemically 



■- IN IN 



I I 



I I I 



I I I I 



CD > 

.It 



z 

■a 

c 

3 

o 
a 
E 
o 
u 



i (T3 

gs -£ 

0 _j 

c O aj 

S o ° 

2 « S 

o Si S, 

Q >o 9- 

1 ™- ° 

O fN -D 



q 



q 



o o o o 



o> a) u) 



£ 
a 
a 



<"0 VO VD UO 

rM ro vo ro co 
co uo vd co en 



vD vD vD 



rM i- i — LT) 

rM m ro rsi 



CTi CO 



N 

E 



ro Ln vD CO 



co o\ ro ro 
o> m co cm 



LO OJ t— 



rM ro ro rM 



5 


CO 


o 


o 




o 


CM 

ro 


CO 


o 

fM 




uri 


oS 


5 


ro 


ro 
ro 



S Ca CO CO 

co co co ro uo uo 

a\ a\ 

o> d CO 

rM t— 



Q\ 0^ CO vO 



ro Oi ro ro 



CO vO 



C 

o 
u 



o 
z 



co a* o «- 



ro ro ro ro 



rM r-i cm 



£ S. 



<T3 

Q- IV 

Q. ro 

03 

h.2 
°- 

(u O 
Q- "D 



PLOS ONE | www.plosone.org 



6 



July 2014 | Volume 9 | Issue 7 | e101432 



Uncovering the Molecular Mechanisms of DB in Colitis Systemically 






R., R2 R3 

25 OH H CH=CH2 

M5 OH OH CH=CH2 

14 h h 0-C2H5 




O 



R1 R2 R3 

4 OH OCH3 OH 

5 OH OH H 
26 OH H h 
M3 H OH OH 
M4 OH OH OH 




R 



R 

37 H 
M8 OH 




R1 R2 R3 R4 R5 

3 H OH OH H OH 

28 H OH OCH3 H OCH3 

M1 OH OH OH H OH 

M2 H OH OH OH OH 

M7 OH OH OCH3 H OCH3 




OH 

41 




R1 R2 R3 R4 R5 

6 OH H OH OH H 
24 OCH3 H OCH3 OH h 
M60CH3 OCH3 OCH3 OH Q H 
27 OCH3 OCH3 OCH3 OH H 



Figure 3. The structures of the 20 potentially absorbed DB constituents and 8 metabolites. 3. 3,5,4'-trihydroxystilbene; 4. 4'-methoxy- 
3'7-dihydroxyflavan; 5. 7,4'-dihydroxyflavan; 6. 2,4,4'-trihydroxydihydrochalcone; 7. norisoboldine; 10. dracaenin A; 13. 3,4-o-dimethylcedrusin; 14. 
ethyl 4-hydroxybenzoate; 16. cochinchinenin; 18. (2S)-7-hydroxyflavanone; 19. 4,4'-dihydroxydiphenylmethane; 21. homoisosocotrin-4-ol; 24. 
loureirin A; 25. 3,4-dihydroxy-allylbenzene; 26. 7-hydroxyflavanone; 27. loureirin B; 28. pterostilbene; 32. dracaenol B; 34. cochinchinenin E; 37. 
diosgenin; 41. bincatriol. M1 and M2 were metabolites of compound 3; M3, M4, M5, M6, M7, and M8 were metabolites of compounds 4, 5, 25, 27, 
28, and 37, respectively. 
doi:1 0.1 371 /journal.pone.01 01 432.g003 



7.3 Pathway enrichment analysis for the major putative 
DB targets and the major known therapeutic targets for 
colitis. We performed pathway enrichment analysis using 
pathway data obtained from the FTP service of KEGG [43] 
(Kyoto Encyclopedia of Genes and Genomes, http:/ /www. 
genome.jp/kegg/, Last updated: Oct 16, 2012). 



Results and Discussion 

1. Optimization of UPLC-ESI-MS/MS conditions 

To obtain reliable chromatographic results and appropriate 
ionization, four mobile phase systems of acetonitrile-water, 
methanol-water, acetonitrile-acid aqueous solution, and metha- 
nol-acid aqueous solution were tested and compared. The results 
suggested that acetonitrile-formic acid aqueous solution was 
superior to the others. Different concentration of formic acid 
aqueous solution (0.05%, 0.1%, and 0.5%) were investigated, 
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Figure 4. The MS/MS spectrum and possible fragmentation pathways of 7,4'-dihydroxyflavan. 

doi:10.1371/journal.pone.0101432.g004 



indicating that use of 0.1% formic acid aqueous solution produced 
the optimal peak shape and reduced peak tailing, as well as 
increasing ion response for most compounds. Meanwhile, the 
process of gradient elution was optimized to obtain better 
separation over shorter time periods. The MS parameters were 
optimized to improve ion intensity. These optimized parameters 
were described in section 2.3. Base peak chromatograms of DB in 
positive and negative ion modes are shown in Figure 2. 

2. Screening and identification of DB using UPLC-ESI-MS/ 
MS 

Recently, the UPLC-Q;Exactive system was reported to provide 
a rapid, reliable and accurate technique for identification of herbal 
constitutions, because of its separation efficiency, high resolution 



and high mass accuracy [44-46] . In the present study, a total of 48 
components of DB were identified. These were mainly flavonoids, 
such as flavones, flavonols, flavanones, chalcones, and isoflavones, 
as listed in Table 1. For some compounds (7,4'-dihydroxyflavan, 
loureirin A, loureirin B, pterostilbene, and dracaenin A), purified 
standards were available and these were identified by comparing 
sample retention time and accurate mass with that of the standard. 
For compounds where standards were unavailable, identification 
was based on accurate mass and tandem mass spectra. The 
molecular formula was established by high-accuracy quasi- 
molecular ion such as [M-H]", [M+CH3COO]", [M+H] + , and 
[M+Na] + within a mass error of 5.0 ppm and fractional isotope 
abundance. The most rational molecular formula was identified 
using chemical or herbal databases, such as Chemspider (www. 
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Figure 5. DB chemical component-putative target network. 

Blue squares refer to the DB compounds and metabolites; yellow 
spherical nodes refer to the predicted targets. 
doi:10.1371/journal.pone.0101432.g005 

chemspider.com) and TCM Database@Taiwan (http://tcm.cmu. 
edu.tw). When several isomers were matched, a compound that 
had previously been identified in DB was considered more likely to 
be correct. Finally, ion fragments were used to provide further 
confirmation of the chemical structure. The structures of the main 
constitutions of DB are shown in Figure 3. For example, 7,4'- 
dihydroxyflavan (compound 5) was firstly identified by comparing 
its retention time and accurate mass with that of the appropriate 
standard. Then, the MS/MS spectrum and possible fragmentation 
pathways of 7,4'-dihydroxyflavan were depicted in Figure 4. 
Only [M+FfJ + and [M-H]~ were observed in positive and 
negative ion modes, respectively. The loss of 128 Da and 
1 36 Da were the characteristic neutral loss of 7,4'-dihydroxyflavan 
in positive and negative ion modes. 

3. Prediction of absorbed DB constituents and their 
metabolites by in silico ADME models 

Quantitative structure-permeability relationships (QSPerR) 
have been established for permeability across Caco-2 monolayers 
(Papp) by the application of new molecular descriptors [47]. 
Among them, the passive Caco-2 absorption model provides a 
useful drug discovery tool by predicting human oral absorption of 
new chemical entities [48-49]. Yazdanian el al. [50] reported that 
compounds with Papp values less than 0.4x10 cm/s exhibited 
very poor oral absorption, whereas compounds with Papp values 
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Figure 6. Putative DB targets-known colitis therapeutic targets protein-protein interaction (PPI) network. (A) The network between all 
targets and other human proteins. (B) The network of hub proteins in network (A). (C) The network of the major putative DB targets and the major 
known colitis therapeutic targets in network (B). Yellow spherical nodes indicate the putative targets; pink spherical nodes indicate the known 
therapeutic targets; purple spherical nodes indicate other human proteins that interact with putative targets or known therapeutic targets. Red 
edges in (C) indicate the PPIs of targets involved in the NOD-like receptor signaling pathway. 
doi:10.1371/journal.pone.0101432.g006 



greater than 7xlO b cm/s had excellent oral absorption. The 
present study employed an in silico passive Caco-2 absorption 
model within the ACD/Percepta software to identify DB 
constituents with excellent oral absorption values (greater than 
7x10 6 cm/s) and the results were shown in Table 1. Subse- 
quently, the ACD/Percepta software P450 regioselectivity module 
was applied to predict the metabolites of these well-absorbed DB 
constituents. These results were listed in Table 2 and the 
structures of the metabolites were shown in Figure 3 . 

4. Analysis of DB pharmacological mechanisms in colitis 

In order to analyze the synergistic effects and pharmacological 
mechanisms of DB in colitis, we constructed a chemical 
component-putative target network, and a putative DB targets- 
known colitis therapeutic targets PPI network. Detailed informa- 
tion relating to these similar drugs and putative targets is provided 
in supplementary Table S3. The performance of this prediction 
method was evaluated in our previous study [34] . 

4.1. DB chemical component-putative target 
network. As shown in Figure 5, this network consisted of 67 
nodes (15 absorbed components, 4 metabolites and 48 putative 
targets) and 120 edges. The mean number of putative targets per 
chemical component/metabolite was 2.53. The absorbed DB 
component 14 and the metabolite M5 of compound 25 (3,4- 
dihydroxy-allylbenzene) both had the highest degree distributions, 
and hit 23 and 14 putative targets, respectively. Since chemical 
components or metabolites with a higher degree in the network 
have been demonstrated to be more pharmacologically important, 
and our data indicated that components 25, 32, and the M5 
metabolite had the shortest paths between each other, they may 
play major roles in this pharmacological network, with synergistic 
effects. 

4.2. Putative DB targets-known colitis therapeutic targets 
PPI network. Putative DB targets-known colitis therapeutic 
targets PPI network was constructed using PPIs between putative 
DB targets and known colitis therapeutic targets. As shown in 
Figure 6A, this network consisted of 2397 nodes and 4005 edges. 
According to the previous study of Li et al. [51], if the degree of a 
node is more than 2 fold of the median degree of all nodes in a 
network, such node is believed to play a critical role in the network 
structure, and can be treated as a hub node. Thus, we identified 
463 hub proteins. Then, the PPI network of these hub proteins was 
constructed using the direct PPIs between them. As a result, this 
network consisted of 463 nodes and 1574 edges (Figure 6B). In 



addition, we calculated four topological feature values of all 463 
nodes in the PPI network of hub proteins according to the 
definition described in section 2.7.2. We defined a hub protein, the 
four topological feature values of which were all larger than the 
corresponding median values, as a major hub protein. After the 
calculation, the median values of the four topological features were 
3, 0.06, 31.30, and 3 for "degree," "betweenness," "closeness," 
and "k value," respectively. Therefore, we determined that hub 
proteins with "degree" >3, "betweenness" >0.06, "closeness" > 
31.30, and "k value" >3 were major hubs. As a result, 106 
proteins were identified as major hubs, including 18 putative DB 
targets and 16 known therapeutic targets for colitis. Detailed 
information on the topological features of these major hubs is 
provided in supplementary Table S4. 

To investigate the pharmacological mechanisms of DB in colitis 
further, we analyzed the direct PPIs between the major putative 
DB targets and the major known therapeutic targets for colitis. As 
shown in Figure 6C, this network contained 34 nodes and 42 
edges. According to the results of enrichment analysis based on the 
KEGG pathway [43], these major target proteins were frequently 
involved in the NOD-like receptor signaling pathway (KEGG ID: 
map04621). The intracellular NOD-like receptor (NLR) family 
contains more than 20 members in mammals and plays a pivotal 
role in the recognition of intracellular ligands. NODI and NOD2, 
two prototypic NLRs, sense the cytosolic presence of bacterial 
peptidoglycan fragments that have escaped from endosomal 
compartments, driving the activation of NF-kappa-B and MAPK, 
cytokine production, and apoptosis [52]. Accumulating studies 
have demonstrated that several members of the NOD-like 
receptor signaling pathway, such as heat shock protein 90 
(Hsp90) [53], interleukin-6 (IL-6) [54], tumor necrosis factor 
alpha (TNF-ot) [55], NF-kappa-B inhibitor kinase alpha (IKKA) 
[56], IKKB [56], and adrenergic receptors [57], maybe associated 
with the pathogenesis of colitis. Of these proteins, Hsp90, AD RBI, 
and ADRB2 were identified as putative DB targets. Moreover, a 
number of putative DB targets also interacted with members of the 
NOD-like receptor signaling pathway, as shown in Fig. 6C. 
These findings suggested that the therapeutic effects of DB on 
colitis may involve this pathway. 

Conclusion 

High-throughput analysis, in silico ADME prediction, and 
network pharmacology techniques have recently emerged as 
powerful tools to provide new insights into the active compounds 
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and the molecular mechanisms involved in the therapeutic actions 
of traditional Chinese medicines. In this study, a rapid, sensitive, 
and high-throughput analytical method using UPLC-ESI-MS/MS 
was firstly developed to characterize 48 constituents and 
accurately identify 24 constituents of DB tablet. Then, 22 
components of DB tablet were predicted to be absorbed by an 
in silico passive absorption model, based on the Caco-2 cell 
monolayer. Eight metabolites were predicted using a computa- 
tional module of P450 regioselectivity. Finally, these compounds 
were predicted to interact with 26 putative targets. Using this 
information, pharmacological networks were built to visualize the 
synergistic interactions between DB components and their targets. 
Furthermore, a PPI network of putative DB targets and known 
therapeutic targets for colitis was constructed to pinpoint the key 
targets and pathways. This approach offered a valuable opportu- 
nity to deepen understanding of the pharmacological mechanisms 
of DB in colitis. 

In summary, this research provided a more accurate analysis of 
DB constituents, compared to using TCM databases. Consider- 
ation of ADME using the passive absorption and P450 metabolism 
modules in silico predicted the in vivo situation more closely. 
Finally, network pharmacology analysis was applied effectively and 
helped to interpret the essence of "synergy" in Chinese medicine. 
This provided a new way to identify the key candidate targets and 
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